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Abstract. In this note we give a numerical expression for the bandwidth bw(P^) of the d- 
product of a path with n edges, P„. We prove that this bandwidth is given by the sum of 
certain multinomial coefficients. We also show that bw(P^) is bounded above and below by 
the largest coefficient in the expansion of (1 + x + ■ ■ ■ + x n ) , with k £ {d, d + 1}. Moreover, 
we compare the asymptotic behavior of bw(P^) with the bandwidth of the labeling obtained 
by ordering the vertices of P„ in lexicographic order. 



1. Introduction 



Given a graph G = (V,E), a labeling f of G is a bijection, / : V — >■ [n], where n = \V\. The 
bandwidth bw{f) of / is defined as 

bw(f) = max |/(u) - f(v)\. 

(u,v)GE 

The bandwidth of G, denoted by bw(G), is the minimum bandwidth over all possible labelings 
of G, that is, 

bw(G) = mmbw(f). 

The bandwidth problem consists of finding the bandwidth of a given graph. Lai and Williams |9j 
provided an end-of-the-century survey containing some of the main results known at the time. 

In this note we give a numerical expression for the bandwidth of the (i-product of a path with 
n edges, P~. The vertices of P% are labeled with elements of the set {0, 1, . . . ,n} rf , and two 
vertices u and v are connected by and edge if and only if u and v differ only at one component 
in which they have absolute difference one (see Figure [2] for an illustration of P^)- 

One of the classical results on the field is the following explicit formula for the bandwidth of 
the hypercube Qd = Pf of dimension d > 0: 

(1.1) bw(Q d ) = fV,!,) 

t=0 VL 2 J/ 

This result first appeared, without proof, in [6j (cf. [71 Corollary 4.4]). The first appearance 
in the literature of a proof for (jl.ip was given by and utilizes the recursive nature of the 
Hales order, which we describe in the Section [2j 

In his pioneering paper [6j, Harper shows that the labeling given by the Hales order, defined 
in the following section, solves the bandwidth problem for Q n . His result was later generalized 
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by Moghadam [10], who showed that the Hales order, under the natural modification of the 
weight function, also solves the bandwidth problem for graphs that are product of paths of any 
length. 

In this note we provide a formula for bw(P^) in terms of sums of coefficients in the expansion 
of (1 + x + x 2 + • • • + x n ) d . In particular, the central multinomial coefficients provide a lower and 
upper bound for bw(P^). Furthermore, we show that bw(P^) is o(bw(L^)) as d — > oo, where 
is the lexicographic order of the vertices of P d . This asymptotic comparison was our original 
motivation to study the bandwidth of the product of graphs. The lex order was proposed to 
us by a colleague as a labeling to approximate the bandwidth for the discrete Laplacian. It 
turns out that the lex order is much worse than the Hales order. This asymptotic comparison 
is carried out in Section 13.31 

Theorem 12.31 and its proof are extensions of \U\ Theorem 2] and its proof. 

For the rest of the paper, n will denote an arbitrary fixed integer larger than one. 



The Hales order lists the elements of V(P~) as follows: u < v if and only if (i) w(u) < w(v) 
or (ii) if w(u) = w(v) and u is greater than or equal to v in lexicographic order relative to the 
right-to-left order of the coordinates. Here, w{x) denotes the weight of x G V(Qd) defined as 
the sum of the components of x. Hales order is often referred to as graded reverse lexicographic 
order. 

For instance, the vertices of P 2 2 = {0, 1, 2} 2 are listed as follows: 00 < 01 < 10 < 02 < 11 < 
20 < 12 < 21 < 22. The order can be obtained by sweeping the hyperplane y + x = from 
(0, 0) to (2, 2), as depicted in Figure [2j Each time the hyperplane intersects vertices of P% , they 
are listed according to the direction indicated in the figure. 



2. Hales order on p£ 




Figure 1. Hales order on P 2 2 = {0, 1,2} 2 



2.1. Recurrence nature of Hales order. The Hales order lists the elements by weight and 
then by reverse lexicographic order from right-to-left (so 2 < 1 < 0). The Hales order exhibits 
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a recursive definition, which we now describe. For < k < nd define the (n + l) d x d matrix 
H% with entries in {0, 1, . . . , n} as follows: 



(2.1) 



at 



■ An,d)- 
^(n,<2) 



A 



(n,d) 



A 



nd ■ 



where the Af 1 '^ are defined below. Here k denotes the column vector whose components are 
all k, with < k < nd, and we let A^' 1 ^ = [i] for < i < n. For d > 2, define 



1) A, 



(n,d) 



[0 



o, 



2) for 1 < k < n - 1, 



A 



(n,d) 



' A (n,d-1) 
^(n,d-l) 



A (n,d-1) 



k 

k- 1 





3) for n < k < nd — n, 



A 



(n,d) 



■ ,(n,d-l) 
A k~n 

An,d-1) 
^k-n+l 



n 



n- 1 



A {n,d-l) q 



4) for nd — n + 1 < k < nd — 1, 



A (n,d) _ 



■ An,d-l) 
A k-n 

An,d-1) 
^k-n+1 



n 



n 



A^/'V k + n-nd 



, and 



5) A 



(n,d) _ 



n n 



n] . 



For instance, 
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[0 0] , A 



1 



(2,3) 



1 
1 , A$ 
10 



(2,3) 






1 



1 

2 



2 

1 1 

1 

2 

1 




A 



3 



(2,3) 



1 2 

1 2 

2 1 

1 1 1 

2 1 

1 2 

2 1 




2 2 

1 1 2 

2 2 

1 2 1 

2 11 
2 2 




12 2 

2 1 2 , and A$ 
2 2 1 



[2 2 2] . 



Notice that the number of rows in A k is the number of terms of degree k in the expansion 
of (l + x + --- + x n ) d . 

The following theorem relates the order of the rows of H%, from top to bottom, to the Hales 
order on P d . 

Theorem 2.1. The rows of H d listed from top to bottom orders the vertices of in increasing 
Hales order. 

Proof We proceed by induction on d. As in the proof of |11| Theorem 1], it is enough to show 

that the rows of A^ 1 '^ are all distinct, have weight k and are sorted, from top to bottom, in 
decreasing lexicographic order from right-to-left. 

The case d = 1 follows trivially. Let us assume that the result holds for all integers at most 
d — 1. Notice that all the rows of are different from each other, since, by induction 

hypothesis, the rows of H^ -1 were. Furthermore, the weight of the rows of 4 M) is k since the 

rows of Am' d ^ have weight m. Finally, notice that the vectors of A^ 1 are listed, from top 
to bottom, in decreasing lexicographic order from right-to-left, as they share the leftmost d— 1 
components with vectors in H^ -1 . □ 

By abuse of notation, we denote the labeling given by the Hales order on P d by H d . In 
particular, 11$ (11) = 5. The Hales order gives the bandwidth of P„, that is, 

Proposition 2.2 (Corollary 3.3, [10J). For d>0, bw{P%) = bw{H%). 

The matrix allows us to find bw(P„) by counting the maximum number of rows between 
two vertices of P d . So we let R(A^ l ' d ') denote the number of rows of A^ 1 "'^. In particular, 
notice that R{A^ L,d ^ > ) is given by the coefficient of the term of degree k in the expansion (1 + 
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x + x 1 H h x n f. That is, 



(2.2) ( d 



ko, k±, . . . , k n 



where the sum is over ko,k±, . . . ,k n satisfying Yl^i = d,^,^ = k. Let R n ^d be the array 
containing the nd + 1 coefficients of (1 + x + x 2 + • • • + x n ) d sorted in decreasing order, that is, 
RnJl] > R n M] > ■ ■ ■> R n Jnd + 1] and let 



n 

K(n,d) = y^R n! d[i], 



that is, 1Z(n, d) is the sum of the n largest coefficients of (1 + x + • • • + x n ) d . 

We now present the main result of this note. 
Theorem 2.3. For d>l, 

d-i 

bw(H d ) = ^2TZ(n,i). 



Proof. Let (u,v) G E(P d ) and assume w(v) > w(u) so w(v) — w(u) = 1, that is, u G Aj^ ,d ' and 
u G ^i+i^ f° r some /c. Then either (i) u and v do not share the d th coordinate or (ii) they do. 

Case (i) turns out not to be relevant in the optimization problem. Let u G [^"^ ^ b] and 



v G [^Jfr b + 1], < b < n. If xi-m denotes the first m components of a vector, then 

j(n,d-l) 
i k-b+l 



ux-i-i = Thus there exists v' G [AfeL? b] satisfying H d (v) - H d {u) < H d (y') - H d (u) 



and (u, ?/) G P d . 



So let us assume that u G [A^!^ ^ b] and w G [A^^, ^ b], < b < n. By Theorem 12. 1} we 
have 

- = J ff^ 1 K (i _ 1 ) - ^(^d-i) + K(u,v,b), 

where R d (u,v,b) is the number of rows between u and v in ff^ that do not have b as last 

,{n,d- 
l k-h 



component. Note that each row counted by R d (u, v,b) is in a matrix [At^i! ^ h] with h G 
{0, 1, . . . , n} \ {6}. There are n such matrices, and since the number of rows of A^ 1 ^ ^ is 



[x k - h ](l + x + • • • + x n ) d -\ it follows that 

bw(H d ) < bwiHt 1 ) + K{n, d-1). 

To verify that the inequality is tight, set b = [n/2\ , c = \n{d — l)/2j and let it G [j4c' d 1 b] 
and v G [A n c ^ 1 b], with iT* -1 ^^-!) - H d ~ x {u x ._ d ^i) = bwiH^' 1 ). Notice that the number of 
rows in H d between u and v that have b as last component is bw(H d ~ l ). Furthermore, by choice 
of c, the rows counted by Rn(u, v,b) correspond to the first n elements of the array R^d-ii f° r 
they correspond to the "middle" n coefficients of the expansion of (1 + x + • • • + x n ) d ~ l (see 
Remark [320, that is, R*(u,v,b) = K(n,d- 1). □ 
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Remark 2.4. If n = 1, then 1Z(n,i) = ^ijj anc ^ one recovers the result of Harper [TJ Corollary 

4.4]. Furthermore, bw(H^) = TZ(n,Q) + lZ(n, 1) = 1 + n, and one obtains a special case of 
Chvatalova's Theorem 1.1] stating that bw(P^ x P^) = min{n + l,m + 1}. 

2.2. The case n = 2. One can compute the bandwidth of P2 utilizing the following lemma, 
whose proof follows directly from (12,2ft . 

Lemma 2.5. fl^) = E^J [ A . k J^- 

Table J] exhibits the value of bw(P*) for < n < 8 and 1 < d < 11. 



Table 1. for 1 < d < 11, 1 < n < 8 



d/n 


1 


2 


3 


4 


5 


6 


7 


8 


1 


1 


1 


1 


1 


1 


1 


1 


1 


2 


2 


3 


4 


5 


6 


7 


8 


9 


3 


3 


8 


14 


21 


30 


40 


52 


65 


4 


6 


21 


48 


91 


155 


243 


360 


509 


5 


12 


56 


172 


404 


831 


1514 


2574 


4085 


6 


22 


152 


617 


1835 


4512 


9655 


18716 


33551 


7 


42 


419 


2289 


8464 


25098 


62474 


138816 


279441 


8 


77 


1169 


8463 


39489 


140059 


408667 


1035692 


2352135 


9 


147 


3292 


32011 


185814 


793765 


2695090 


7823236 


19956152 


10 


273 


9338 


120439 


880174 


4499506 


17887694 


59241709 


170376339 


11 


525 


26641 


460813 


4191494 


25788102 


119335481 


452484637 


1461956288 



3. Bounds involving multinomial coefficients 



In this section we show that bw(P^) is bounded below and above by multinomial coefficients. 
We first discuss some combinatorial background. 

3.1. Extended Pascal's triangles. Let C n (d, k) denote the coefficient of x k in the expansion 
of {l + x + --- + x n ) d . 

The array 

C n (0,0) 

C n (l,0) C n (l,l) 

C n (2,0) CW(2,1) 

C n (d,Q) C n (d,l) ••• C n {d,dn). 

is called an extended Pascal's triangle (or generalized Pascal's triangle), and we denote it by 
T d . That is, T d is the array where the element in the (£ + l) st row and (k + l) st column is 



C n (l,n) 
C n (2,2n) 
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C n (£, k), with < I < d and < k < in. These extended Pascal's triangles satisfy properties 
akin to those of Pascal's triangle. We list the properties that will be utilized in the sequel. 

Lemma 3.1. Let d,n,k be integers with < k < d. Then, 

n 

(i) C n (d, k) = Y^ C n (d —l,k — j), where C n (a, b) = 0ifb<0. (See gj Equation (2)].) 

j=0 

(ii) C n (d, k) = C n (d,nd — k), that is, each row ofT d is symmetric. (See [H Equation (4)].) 

(hi) The sequence {C n (d, fc)}Jjf. is log-concave. In particular, each row of T d is unimodal. 
(See [21 Theorem 7, Corollary 8].) 



We refer to Lemma l3.1( i) as the extended Pascal's identity. Some remarks are in order. 



Remark 3.2. Combining the second and third part of the lemma, we conclude that the sequence 
{C n (d, k)}f =0 is unimodal and has a peak at k = [(nd + l)/2j ; furthermore, just as in the case 
n = 1, this sequence can have at most two peaks. The lemma also implies that the largest 
i elements of each row of T d are the "middle" i elements in each row, that is, the elements 
located in positions j, where 

' [{nd + 1)/2J - [i/2\ <j< [(nd + 1)/2J + [i/2\ 
[(nd+ 1)/2J - i/2 <j< [(nd+ 1)/2J + i/2 - 1 
[{nd+ 1)/2J - L*/2J < J < L("d+ 1)/2J + LV 2 J 
s L(nd + 1)/2J - i/2 <j< [(nd+ 1)/2J + i/2 - 1 



if d is even and i is odd, 
if c? is even and i is even, 
if d is odd and i is odd, and 
if d is odd and i is even. 



3.2. The bounds. Using the notation of Section [2.11 we write R(Aj^' ) instead of C n (d,k). 
We furthermore denote the highest coefficient in the expansion of (l + x + --- + x n ) d by M(n,d). 
From Remark l3.2[ the first i elements of Rnd are those that correspond to the "middle" i degrees 
of the expansion of (1 + x + • • • + x n ) d . Therefore, the extended Pascal's identity yields 

n+l 

(3.1) M(n, d) = Y, Rn,d-M = d-l) + R n ,d-l[n + 1]. 

i=i 



Now we prove that bw(P d ) is bounded below and above by central multinomial coefficients. 
The key element in the proof is Equation (|3.ip . 



Lemma 3.3. For n,d>0, M(n, d) < bw{P d ) < M(n, d+l). 



Proof. We prove both inequalities by induction. The base case d = 1 is easily verified to hold. 
Let us suppose that M(n,d- 1) < bw(P d ~ 1 ) < M(n,d). Then, 

M(n, d) = K(n, d - 1) + R^d-x [n + 1] 

< n(n,d-l) + M(n,d-l) 

< Kin^d-^ + bwiPt 1 ) 
= bw(P d ). 
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For the upper bound, notice that 

bw(P d ) = bwiPt 1 ) + d-1) 
< M(n,d) +n(n,d- 1) 

n 

= M(n,d) + J2 R n,d-i[i\ 
1=1 

n+1 

<M{n,d) + J2 R nM 

i=2 

= M(n,d + 1), 

as desired. We used the fact that YH=i ^n,rf-lH — X^i=2 -^n,d[*]) f° r the elements of i? ra ^ are 
sums of elements of R n d-i by the extended Pascal's identity. □ 



3.3. Comparison with the lex order. Our interest in the bandwidth problem started when 
a colleague, John Hubbard, asked whether the lexicographic order provided an asymptotically 
good approximation to the bandwidth for the discrete Laplacian. We point out that one can 
obtain information regarding the bandwidth of a graph via the discrete Laplacian. For example 
the Laplacian spectrum of a graph G can be used to give a lower bound on the bandwidth of 
G (see, e.g., [HH]). In this section, we compare the asymptotic behavior of the lexicographic 
order (from left-to-right) and the Hales order. 

Let denote the label given by ordering the vertices of P d lexicographically (from left-to- 
right). Then it can be proved inductively that 

bw{L d n ) = {n + l) d -\ 

For the rest of the section, we consider n to be fixed. We will show the following proposition. 
Proposition 3.4. bw(H d ) is oQrw^L^)) as d —> oo. 



Proof. We can think of independent random variables uniformly distributed 

on {0, 1, ... , n}. The expected value E(x l ) of each x l , < i < n, is ^ and so 

the variance o~ 2 (x l ) of each x l is 

, 2 n 2 n 2 + 2n 



1 

E((z*) 2 ) - (E(z*)) 2 = -— £ 



n + 4 12 

i=0 



Hence the mean and variance of the (d+ l)-fold convolution are n(d + l)/2 and (d + l)(n 2 
2n)/12, respectively. By the central limit theorem, the coefficients of 



p(d) :- 



n \ d+l 

Xq + ■ ■ ■ + X 

n+1 



are (asymptotically) normally distributed. Furthermore, since the coefficients in the expansion 
of p(d) are log-concave, O Lemma 2] gives that these coefficients satisfy a local central limit 
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theorem around the coefficient M(n, d + l)/(n + 1) . That is, 

M(n,d + 1) 1 _j 

as \j \p2/n is the largest value of the density function of the standard normal distribution (hence 
we obtain an overestimate). Therefore, 



M(n,d+1) ~ (n + l) d+1 , 

Therefore, 
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vr(d+ l)(n 2 + 2n)' 



,. ^ .. M(n,d+1) , 

lim ; " y < hm ~ (n + l) z 



d^oobw{L d n ) ~ d^oo (n + iy- 1 K ' V vr(d+ l)(n 2 + 2n) 
which tends to zero as d — )• oo. □ 



4. Further directions 



It would be interesting to obtain some kind of bound for the bandwidth of the product of 
paths of different length. Moghadam [TU] showed that the Hales order is optimal for these 
graphs as well, so finding a recursive procedure to gives the Hales order for these graphs might 
prove helpful to obtain bounds. 

Finally, it would also be interesting to obtain an explicit formula for the antibandwidth of 
the product of paths. The antibandwidth of a graph G = (V,E) is given by 

abw(G) = max min \f(u) — f(v)\. 

f u.vGE 

That is, abw(G) maximizes the minimum absolute value of the difference of labels of any 
adjacent pair of vertices, over all labelings of G. In the authors provide an explicit formula 
for abw(Pf), but to the best of our knowledge no such formula exists for the product of paths 
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